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Search for Lepton Flavor- Violating r — > fij Decay 

K.Inami a , T.Hokuue a , T.Ohshima a for the Belle collaboration 

a Department of Physics, Nagoya University, 
Furo-cho, Chikusa-ku, Nagoya, 464-8602, Japan 

We search for the lepton flavor- violating r — > /ry decay using 29.7 million r pairs accumulated by the Belle 
experiment. The main background sources are found to be r pairs with r — » /xz/i/ decay and radiative /i/i events. 
One event is observed in the signal region, while 2.5 ± 0.6 background events are expected. A preliminary upper 
limit Br(r — > fAj) < 6 x 10~ 7 at the 90% confidence limit is obtained. 



1. Introduction 

Charged lepton flavor-violating decays, such as 
t — » fj/y, ey, /i/i/i and (J, — > ej, are forbidden in 
the Standard Model (or highly suppressed even 
if we consider neutrino mixing). However, new 
physics beyond the Standard Model allows lepton 
flavor violation (LFV) decays. In supersymmet- 
ric models, left-right symmetric models and oth- 
ers |lj , the branching ratio of the r — > /17 decay, 
Br(r — » /17), is predicted to be 10 5 " 6 times higher 
than Br(/i — > &y) because of the large t mass. 
Also, some SUSY models predict Br - O(10 -7 ), 
which can be reachable by B-factory experiments 
of high luminosity. 

The current upper limit of Br(T — * ^7) is 1.1 x 
1CT 6 (90%CL), measured by the CLEO experi- 
ment |Q. Our first analysis with 10 fb _1 of data 
resulted in a slightly better upper limit of 1 x 10 -6 
(90%CL) ||. Also, the BaBar group recently re- 
ported an upper limit of 2.0 x 10~ 6 (90%CL) us- 
ing 56 million t pairs at the ICHEP2002 confer- 
ence [§. 

Data is accumulated with the Belle detector || 
at the KEKB accelerator KEKB is an asym- 
metric e + e~ collider with a center-of-mass en- 
ergy of 10.58 GeV. Its current peak luminosity is 
7.4 x 10 33 cm _2 s _1 and the total integrated lumi- 
nosity amounts to 89.6 fb _1 as of summer, 2002. 
Here we present an analysis using 32.6 fb" 1 of 
data, corresponding to 29.7 million r pairs. 



2. Event selection 

In order to determine the selection criteria, a 
Monte Carlo (MC) simulation is performed us- 
ing the KORALB generator for r-pair sam- 
ples and a QQ generator J^] for other physics pro- 
cesses, as noted below. GEANT |J is used for the 
full detector simulation. We prepare two kinds 
of r pair samples: signal and r-pair background, 
and several kinds of background samples: B°B , 
B + B~, continuum, Bhabha, /i/j, and two-photon 
processes of eeee, ee/x/j,, eeuu/dd, eess and eecc. 
In the signal MC samples, one r decays to fv) and 
the other r decays generically. The angular dis- 
tribution of two-body r — > /17 decay is assumed 
to be uniform in the r's rest frame. 
A candidate event is required to be 
(^7) + [(a charged particle but not fx) 

+ (> O7) + missing], 
which we hereafter denote as "/z n °V" events. We 
therefore select events with "2 charged-tracks and 
>1 7" with the expected accompanying missing 
momentum. Below, kinematical variables are de- 
fined in the laboratory frame, while those defined 
in the center-of-mass frame of the system are 
indicated with an index "CM" . Two oppositely 
charged particles are required to have a transverse 
momentum of pt > 0.1 GeV/c and a momentum 
of p < 4.5 GeV/c to reject Bhabha and /z/x 
background events. For the signal, a muon track 
is required to satisfy —0.819 < cos^ < 0.906 
and p^ > 1.0 GeV/c, and to have a /^-probability 
of more than 0.90 from the KL/muon detector 
(KLM). The signal photon is selected by requir- 
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Figure 1. Comparison of some kinematical distributions among data (dots), r pair background MC (open 
histogram) and signal MC (hatched histogram). The arrows indicate the selection criteria. 



ing -0.866 < cos6> 7 < 0.956 and E 1 > 0.5 GeV. 
The charged particle on the tag side is required to 
not be a muon within —0.819 < cos# t ag < 0.906. 
In order to remove the background from radia- 
tive n/i, two-photon and other non-r processes, 
we apply the following criteria. The opening an- 
gle between the two charged particles be greater 
than 90°, and for the missing particle, p m i ss > 0.4 
GeV/c and 17° < 6 miss < 150°. 

To remove the r-pair background, we demand 
the opening angle between the /i and 7 of the 
signal candidate to be 0.4 < cos 8™ < 0.8 
and the opening angle between the charged par- 
ticle in tag side and the missing momentum to 
be cos6>£^ miao > 0.4. The distributions for 
these variables are shown in Fig. [|. A further 
selection is applied for the missing momentum 
(Pmiss) and the missing mass squared (m miss ) as 
Pmiss > -5m^; - 1 and p miss > 1.5ml liss - 1, as 
shown in Fig. |[ With these selections, the sig- 
nal efficiency is evaluated to be 12.8% by MC and 
the remaining rate of background is5.7xl0 _6 for 
r-pairs and 2 x 10 -6 for the radiative \x\x. 

The signal region on the M\ m -AE plane is de- 
termined using the signal MC while assuming a 
uniform background distribution around the sig- 
nal peak, where Mi nv is the invariant mass of the 
/1-7 system and AE is the energy difference be- 
tween the /i-7 system and the beam energy in 



the CM frame. Although we used a box shape to 
define the signal region in previous analysis, we 
now employ an elliptic shape in order to enhance 
the sensitivity, because the signal distribution has 
a tail due to the initial-state radiation and pho- 
ton energy leakage from the Csl photon detector. 
Fig. |^ shows an ellipse with the highest signal-to- 
background ratio for a total detection efficiency 
of 9.0%. 

With all selections, 69 events remain within the 
regions called "Area", defined to be 1.5 < M; nv < 
2.0 GeV/c 2 and -0.5 < AE < 0.5 GeV. One 
event is observed in the signal region, as plotted 
in Fig. |. 

3. Background estimate 

It is not straightforward to use a sideband to es- 
timate the background in the signal region, since 
it largely depends on the sideband allocation for 
small samples, especially in the case of a struc- 
tural distribution. 

The backgrounds from r pairs distribute in the 
AE < region, as shown in Fig. |I| and could 
possibly contaminate the signal region. While 
a negligibly small portion of r pairs is found at 
AE > 0, the data extend to < AE < 0.2. Our 
study finds these to be radiative events, one in 
which muon is not identified due to the KLM in- 
efficiency Fig. ^(b) shows the events that passed 
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Figure 2. Correlations between p m i ss and mf niss for (a) the signal MC and (b) the r-pair background MC. 
The lines show the selection boundaries applied in this analysis. 
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Figure 3. Mi nv -vs-AE plot for events which 
passed all selections. The data events are indi- 
cated by dots and the MC events by boxes. The 
ellipse shows the signal region. 



the same selection criteria as the signal, with the 
tagged track required to be a muon. We call these 
"/i/i" events. The muon identification inefficiency 
is estimated to be about 8%. 



The background rate is therefore calculated us- 
ing 

^Vbg = A^ n °^(MC) 

+{iV w (data) - TV^(MC)} x [j—^j > M 

where F'^lMC) is the number of "/i not /i" 
events estimated by MC. The background is 
mostly r-pair background and a small portion of 
the continuum. TV w (data) is the number of "/i/j," 
events extracted from data. The last term is in- 
cluded to avoid double counting. The resultant 
background rates of the individual contributions 
in Eq.(jl|) are listed in Table [l]. Only 5 contin- 
uum events (8% in the Area) remain in "/it 110 '//" 
samples for 32.6 fb _1 luminosity in the Area, but 
outside of the signal region. The contribution of 
the fourth term is not significant. It is only 1.9% 
in the Area. 

A comparison between the number of the ob- 
served events and the expected backgrounds is 
given in Fig. The center box is the signal re- 
gion used in the previous analysis. We find a good 
agreement between them. We observe 69 events 
with an expected background of 62.5 ±3.5 events. 

Comparisons of the AE and M- mv distribu- 
tions between data (dots) and the expected back- 
ground (open histogram) are also shown in Fig. 
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Figure 4. M inv -vs-A£; plot of (a) the r pair background MC (140M) and for (b) A^(data) (32.6 fb _1 ). 
The A MAI (data) component contributes to the background through the muon identification inefficiency of 
about 8%. 

Table 1 

Number of the observed events and the expected backgrounds for individual components of Eq. ([I]) for a 
luminosity of 32.6 /ft -1 . 
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Figure 5. Number of observed events (left) and 
the expected backgrounds (right) in each region. 



While a uniform distribution is seen for Mj nv , a 
two-bump structure is found in the AE distri- 
bution. The two bumps corresponds to t pair 



backgrounds and radiative processes. Fig. fj] 
shows the and E 7 distributions of the signal 
candidate for the same samples of Fig. ||. Both 
bumps in the high-p^ and the low-_E 7 regions are 
due to the remaining radiative /i/i samples. As 
can be seen in Figs. |^ and ^, the distributions 
of the expected backgrounds are very consistent 
with the data. 

4. Result 

In the signal region, we find one event, as shown 
in Fig. ||. The number of background events is 
evaluated to be 2.5 ± 0.6. The upper limit is 
obtained using the Bayesian approach (following 
refs [9, Jin] and 111)) with the equation 



e- & °E;;icA7«! 



0.1, 



(2) 



where sq is the upper limit on the signal at 90% 
confidence, 6q is the number of expected back- 
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Figure 6. (a) AE and (b) M- mv distributions for the surviving events in the Area. The points with error 
bars are data, the open histogram is the expected background evaluated from Eq.(|l]), and the hatched 
histogram is the signal MC (Br(r — > (J,j)=2 x 1CT 6 ). 
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Figure 7. (a) p^ M and (b) E^ M distributions of the events that survived all selections. The points with 
error bars are data and the open histogram is the expected background evaluated from Eq. ([!]). 



ground events and no is the number of observed 
events. We obtain s = 4.1. 

The systematic uncertainties of the detection 
efficiency are listed in Table 0. The conservative 
and preliminary evaluation gives a 6.7% error in 
total. 

Finally, we obtain a preliminary result for the 
upper limit of Br(r — > fij) using 29.7 million r 
pairs with e =9.0% detection efficiency, as 

Br(r -> ^7) < 6 x 10~ 7 . 
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Table 2 

Systematic uncertainties in the detection effi- 
ciency. The total uncertainty is evaluated by 
adding them in quadrature. 





uncertainty (%) 


track rec. eff. 


2.0 


photon rec. eff. 
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cut selection 
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MC statistics 
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